Exact solution for eigenfunction statistics at the center-of-band anomaly in the 

Anderson localization model 
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An exact solution is found for the problem of the center-of-band (E = 0) anomaly in the one- 
dimensional (ID) Anderson model of localization. By deriving and solving an equation for the 
generating function $(«,(/>) we obtained an exact expression in quadratures for statistical moments 
I q — (\ipE(r)\ 2q ) of normalized wavefunctions tpE(r) which show violation of one-parameter scaling 
and emergence of an additional length scale at E ~ 0. 
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-Introduction. Anderson localization (AL) enjoys 
an unusual fate of being a subject of advanced re- 
search during a half of century. The seminal paper by 
P.W.Anderson 1 opened up a direction of research on the 
interplay of quantum mechanics and disorder which is of 
fundamental interest up to now 2 . The one-dimensional 
(ID) tight-binding model with diagonal disorder -the An- 
derson model (AM)- which is the simplest and the most 
studied model of this type, became a paradigm of AL: 
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Cj+l + c. 
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In this model the hopping integral is deterministic U = 
t = 1 and the on-site energy £j is a random Gaussian 
variable uncorrelated at different sites and characterized 
by the variance ((fe,) 2 ) = w. 

The best studied is the continuous limit of this model 
in which the lattice constant a — s- at to? remaining 
finite 3 -"^. There was also a great deal of activity 5 aimed at 
a rigorous mathematical description of ID AL. However, 
despite considerable efforts invested, some subtle issues 
concerning ID AM still remain unsolved. One of them 
is the effects of commensurability between the de-Broglie 
wavelength A e (which depends on the energy E) and the 
lattice constant a. 

It was known for quite a while^ that at weak dis- 
order w <C 1 the Lyapunov exponent takes anomalous 
values at the ratio / = equal to \ and -| (compared 
to those at / beyond the window of the size w around 
/ = 2 ano - / = 3)- The Lyapunov exponent sharply 
decreases at / = | (which is usually associated with in- 
creasing the localization length) and may both increase 
or decrease at / = -| depending on the third moment 
of the on-site energy distribution 7 . More recently^ it 
was found that the statistics of conductance in ID AM 
is anomalous at the center of the band that corresponds 
to / = |. We want to stress that all these anomalies 
were observed for the AM Eq.([T]) in which the on-site 
energy is random. This Hamiltonian does not pos- 
sess the chiral symmetry!^ which is behind the statis- 
tical anomalies at the center of the band E — in the 



Lifshitz model described by Eq.([T]) with the determinis- 
tic £i — and a random hopping integral t%. Thus the 
statistical anomaly at / = ^, \ raises a question about 
a hidden symmetry that does not merely reduce to the 
two-sublattice division 2 -^^. 

A similar phenomenon may occur in dynamical sys- 
tems. An elegant analogy between the ID localization 
and the classical system of kicked oscillator was studied in 
Refi^i. According to this analogy the energy-dependent 
de-Broglie wavelength A^ is encoded in the frequency of 
the oscillator and the lattice constant a determines the 
period of the 6- function "kicks" of the external force, 
their amplitude being proportional to disorder. 

The interest to one-dimensional Anderson localiza- 
tion is greatly increased recently after several groups re- 
ported about successful experiments on localization of 
cold atom o 12 i 13 , where even tiny details of localized wave- 
functions were observed. Kicked rotors and kicked oscil- 
lator can also be realized in systems of cold atomsii. 

There are numerous questions concerning physics be- 
hind the anomalies. One of puzzles is the sign of the 
variation of the Lyapunov exponent which corresponds 
to weaker localization at / — ^ . Such a tendency can be 
considered as a remnant of the chiral symmetry spoiled 
by fluctuating on-site energy. There is, however, a com- 
pletely different view on the problem which predicts the 
stronger localization at the band center. It involves the 
notion of Bragg mirrors^ created by disorder realiza- 
tions with alternating on-site energies which double the 
period, at least locally. A possible resolution of this 
conflict between different mechanisms of the center-of- 
band anomaly could be a typical wavefunction sketched 
in Fig.l. It contains two length scales: one of them £\ oc 
which is somewhat larger than the localization length £ 
away from the anomaly, is due to remnants of the chiral 
symmetry, while the other, much smaller one d <C £\ oc 
(but d 3> a) , is due to the formation of the Bragg mirror 
fluctuation. If the weight of the narrow peak p <C 1 is 
small, the statistical moments I q — (\ip(r)\ 2q ) of the nor- 
malized wavefunctions "0( r ) with relatively small q will 
follow the standard^ 7 - behavior (L is the length of the 
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FIG. 1: The cartoon of the standard (blue) and the E = 
typical state (red). At the anomaly the second scale d 
emerges. 



chain) : 

I^ = {q-l)\l\-«L-\ (2) 
while at large q, the higher and more narrow peak will 



dominate in I qi leading to the d - ^ -1 * 1 instead of £7^ ^ 
behavior of moments. The simplest crossover between 
the two regimes is described by 
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Eq. <j3j) can be obtained from the following qualitative ar- 
guments. The average moment is the sum of two con- 
tributions. The first one leading to the standard mo- 
ment Eq.©, is equal to the q-th. power of the typical 
amplitude IVK 1 ")! 2 ~ 1/Aoc inside the localization ra- 
dius but outside of the narrow peak, multiplied by the 
probability ~ £\oc/L that the observation point r falls 
inside this region. The second contribution (p/d) q d/L 
arises when with small probability d/L the observation 
point falls inside the narrow peak where the amplitude 
|^?(r)| ~ p/d. It is this contribution which corresponds 
to the d-dependent term in Eq.Q. 

In general, the information about a typical shape of lo- 
calized wave functions is encoded in statistical moments 
I q . In this Letter we solve exactly the problem of statisti- 
cal moments I q at E — for the ID Anderson disordered 
chain Eq.([T]) of the length L — > oo and show that the 
behavior Eq.Q indeed emerges. 

- Generating function (GF), moments of \ip\ 2 and the 
probability distribution function (PDF) of phase. Mo- 
ments I q of normalized eigenfunctions with integer q > 1 



I,(r) = (|Mr)| 2 *) = 
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dzz q 2 $ r _i(z, 4> ~ k) §N-r(z, ~4> - k), (E = 2cosk). (4) 



can be expressed^ in terms of a generating function rived elementary-i&starting from Eq.([T]) as well as using 
<&j(u, (/)) on the lattice site j. The starting point of our the super-symmetry method^: 
analysis is the recursive equation for GF which can be de- 
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This equation is exact and holds both for weak and strong 
disorder controlled by the parameter w for any energy 
E = E(k) parametrized by E(k) = 2cosfc. 

The way the variables z and <j> enter Eq.Q suggests 
their physical meaning^: they determine the values of a 
wave function tp(i + 1) and on a link {i, i + 1}: 



ip{i) = Vzi cos(^i), tp(i + l) = y/z. 



k). (6) 



It is remarkable that both the "elementary"— and the 
"super-symmetric"— derivations of eigenfunction statis- 
tics involve naturally the two link variables Eq.([B]). In 
contrast to the moments /„, the Lyapunov exponent 
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depends on only one of the two sets of variables, <fi 
(ip(i + l)/ip(i) = cos(0j — k)/ cos(^j)), which determines 
completely its statistics^. That is why the problem of 
moments is more complicated and more general than that 
of the Lyapunov exponent. 

The integrand in Eq.© is bi-linear in $. This effec- 
tively takes into account the boundary conditions at the 
two ends of the chain 1 ^ which is necessary to describe 
the normalized ezgenfunctions. In contrast to that in the 
problem of Lyapunov exponent^ one considers essentially 
a semi-infinite chain and does not require of the solution 
to Eq.© to be an eigenfunction. 

Despite the fact that &j(z, <j>) is not the joint PDF of 
z and 0, its descender $j(z — 0, 4>) is the PDF of phase: 



in ajl we obtain for k = irp/q: 



$,(z = 0,0) =P J (0), 



P,(0)d0 = l. (8) 



This statement can be formally proven 1 ^ but the key 
properties of PDF, the positivity of Pj(<fi) and the con- 
servation of normalization, are easily seen directly from 
Eq.© and the boundary condition $o(0) = S((f) — 7r/2). 

-Evolution equation for weak disorder. Eq. © is valid 
for an arbitrary strength of disorder. However, the 
anomaly we are going to study is sharp only at weak 
disorder and is rounded off as disorder increases. For 
weak disorder when the localization length i = sm 
is large compared to the the lattice constant a one can 
reduce Eq.© to a partial differential equation (DE) 
of the Fokker-Planck type, where the coordinate x = 
ja/£ along ID chain plays a role of time and the two- 
dimensional space of variables u — £z and stands for 
the coordinate space. 

In the first order in a/^ < 1 one obtains by a proper 
expansion in Eq.©: 



$j+i(tt,0) = (l + j [£(0)-ci(0)u]) $i(u,0-fc), 



(9) 



where the evolution operator ^1 + | C(4>) — ci(0) u 
contains the differential part 

t{4>) = c 2 {^)u 2 dl + c 3 (4>){ud u -\) (10) 
+ C4(0) ud u d^ + c 5 (0) 9^ + ce((f)) d\ 

and Ci(cj)) are certain linear combinations of 1, sin(20), 
cos(20), and sin(40), cos(40). 

The formal reason for the center-of-band anomaly at 
k = ? (as well as of the weaker anomaly at any k — irp/q, 
where p, q are positive integers) is the shift by k of the 
argument in r.h.s. of Eq.©. Because of this shift 
and the periodicity $j(u,4>) = $., (it, + 7r) , one has 
to apply the evolution operator q times in order to get 
a closed recursive equation which expresses <&j+ q {u, 0) 
in terms of $>j(u,<j)) and its derivatives. For weak dis- 
order and not very large q <C £/a, one can expand 
&j+q — &j ~ (aq/l) d$(u, 0; x)/dx, where we introduce a 
function $(u, 0; x) = ^i x / a (u, 0) of a continuous dimen- 
sionless coordinate x — ja/t. Thus in the lowest order 
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The sum over s arises because the small corrections to 
the evolution operator proportional to ajl add up in the 
product of q evolution operators, each time entering with 
a shift Ci(<fi) —> Ci(<f>—k) according to Eq.©. The crucial 
point for emergence of anomaly at fc = ^ (q = 2, p = 1) 
is the identity: 
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0, q>2 
v qe 4 ^, q = 2 

(12) 

One observes that at k = irp/q with all q but q = 2 the 
0-dependent terms disappear from the r.h.s. of Eq. ([lip . 
At fc = -|, however, one obtains the anomalous, 0- 
dependent, evolution equation. It appears to have a nice 
SL{2) group structure: 



d x <f>(u, 0; x) = {L\ + L\ - u} $(u, 0; a;), 
L\ = cos9de + sm9ud u , L^ — de, (9 



(13) 



sin 9 dg 



where L\ and L3 and L2 = [L^,L{\ - 
cos9ud u form a closed sl(2) algebra. 

Note that Eq. flT3")) contains all the known particular re- 
sults. For instance, omitting all the 0-dependent terms 
one obtains the standard equation for GF away from the 
anomalyii which allows for the 0- and x-independent 
(zero-mode) solution&i 



(14) 



Alternatively, in agreement with Eq. © , by setting u = 
in Eq. (|13p one arrives at the second order ordinary DE 
for the non-trivial phase-distribution function i"b(0) at 
the k = ^ anomaly with the zero-mode solution: 



P O (0) = $(O,0) = 
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y/3 + cos(40) ' 



c = 
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(15) 



resulting in the anomaly of the Layapunov exponent*^: 



r 2 (?) 

V4/ -0.9139. 



ra = r <1+cosW))PoW - 

(16) 

Derivation of Eq. (ll3|) . and its exact solution is the main 
result of this Letter. 

-Separation of variables and the zero-mode solution. 
The variables u and are entangled in Eq. ([T3|) . How- 
ever, there is a hidden symmetry which allows to sepa- 
rate variables in this equation, provided that the term 
d x &(u, <p]x) — 0. This zero mode solution is sufficient 
to describe anomalous eigenfunction statistics in a very 
long chain L ^ £ far from its ends. 
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' Correct variables" £ and n are suggested by Eq. ([6]) : 
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Denning also the "correct function 



(Z + v) 



*(vfov), <t>(£,n)), 



(17) 



(18) 



one casts the zero-mode variant of Eq. ()13|) in the form of 
the Schroedinger equation: 



H(0+H( V ) * = 0, ff(£) = -# 



16 £ 2 



1 

(19) 



Note that the singular operator H(£) is not Hermitian 
for generic wave function. Its spectrum is continuous 
and, in general, complex. The zero-mode solution cor- 
responds to a zero sum of the two eigenvalues (±A) of 
the ID Hamiltonians H(£) and H(rf). Thus the solu- 
tion to Eq. (|19p emerges as an integral over a continu- 
ous variable A oc 1/ a/A which can be taken real without 
loss of generality The integrand involves the prod- 
uct of two eigenfunctions x I'a(£) and ^_a(^), and an ar- 
bitrary function C(A). Yet, one can find this function 
C(A) uniquely^ using the conditions of (i) smoothness 
of <ft) at <j> = and cj> = ? and (ii) normalization of 
the phase distribution function Pq(4>) — $(u = 0, </>): 
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where W a i(x) is the Whittakcr function, T(x) is the 

Euler Gamma- function, and e = e™/ 4 , e = e~ l7T / 4 . 

Eq. (f2"0)) is the main analytic result of the Letter. 

-Moments of normalized eigenfunctions. A convenient 
way to present the results is to plot the reduced moments 
R q =I q (E = 0)/I q (E^0), 



R„ — C„ 



du 
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0). (21) 



Here I q (E ^ 0) = L" 1 (q - 1)1 l x ~^ arc the moments 
away from the anomaly, where $(it, </>) = $>( st \u) is given 

Using the solution 



7T4" 



by Eq.QJJ), and C q = (g _ 1)!(g _ 2) , ■ 
Eq. (|20|) we evaluated the reduced moments i? g numeri- 
cally up to q — 10. The results are given in Fig. 2. One 
can see that at E = the moments R q « (^/^ioc) 9_1 with 
small q follow Eq.Q, albeit with a localization length £\ oc 
larger than that away from the anomaly. The best ex- 
ponential fit of moments with q < 6 gives l\ oc /l ~ 1.252. 
This reflects the same tendency as Eq. pi)]) . However, 
larger moments are significantly greater than the predic- 
tion of one-parameter scaling Eq.([2]). The excess factor 
S q = I q (E = 0)/I q st ^ which should be compared with 
that in the square brackets of Eq. © , is plotted in the in- 
sert of Fig. 2. A comparison with Eq.Q shows a very sat- 
isfactory (for a crude qualitative interpretation in terms 
of two scales sketched in Fig.l) agreement for moments 
up to q = 10 which can be interpreted as an emergence of 
a very narrow (but still much wider than the lattice con- 
stant) peak in an "average" eigenfunction at the anomaly. 

In conclusion, we solved exactly the problem of sta- 
tistical moments I q of the amplitude IV'.eMI 2 of random 
wave functions in the ID Anderson model at energies 
E w 0. It is shown that the statistics of such wave- 
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FIG. 2: (color line)Reduced moments R q (red points) in the 
log-linear scale. The dashed line is the exponential fit R q = 
(e/eioc)"' 1 with l loc (E = = 1.252. In the insert: The 
excess factor (red points). The dashed line is a fit by Eq.@ 
with p « 0.01 and k oc (E = 0)/d « 160. 



functions is anomalous which anomaly does not merely 
reduce to the variation of the localization length or the 
Lyapunov exponent. The enhancement of the localiza- 
tion length £\ oc /£ ~ 1.252 derived from I q cx il~ c 9 with 
q < 6 is different from that obtained from the inverse 
Lyapunov exponent j(E ^ 0)/j(E = 0) « 1.094. This 
fact together with the anomalous enhancement of mo- 
ments with large q > 6 implies a significant change of the 
form of the typical eigenfunction at E w which requires 
more than one characteristic length for its description. 
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